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Abstract 

Based on the measured transverse mass spectra of ir~ , K~ and p at the RHIC 
energy y / Sj\F/v = 130GeV, di-omega productions from baryon-baryon reactions in 
hadronic matter are studied. Results about the (fif2) + total number per event 
show that the deeply bound state (S7r2) + can be observed at RHIC energies. 
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1. Introduction 



Data measured at RHIC experiments have revealed intriguing new features in ul- 
trarelativistic nuclear collisions [1] . Various theoretical models have undergone challenges 
from the experimental data and turn out to be adjusted. The unexpected large val- 
ues of the elliptic flow signal v 2 can be accounted for by hydrodynamic calculations [2] 
indicating perhaps complete early-time thermalization in central Au+Au collisions at 
v /sjvtv — 130GeV per nucleon pair. This stresses the collective flow strongly affects 
hadrons in the hadronic matter produced at the center-of-mass energy. Therefore, any 
further predictions relying on the hadronic matter must involve available hydrodynamic 
behaviors. 

One of the interesting topics in hadronic physics is the search for dibaryon bound 
states. Since dibaryon is a kind of new matter, it can provide a good place to examine 
the short-range quark-gluon behaviour of QCD theory and also would open a new area 
for studying many new physical phenomena we have not known before. Even though 
the H dibaryon predicted by Jaffe [3] has not been confirmed in hadron-hadron collison 
experiments, recently, the search is extended to Au+Au collisions in the AGS Experiment 
896 [4]. Dibaryons involving A and nucleon were predicted in an adiabatic method [5]. 
Productions of short-lived dibaryons with strangeness, for instance, E°p, S°A, 
have been calculated for Au+Au collisions at y/s NN = 200GeV and may be identified as 
peaks in the invariant mass spectra of final two baryons yielded from dibaryon decays [6] . 
Our present work is based on a dynamical investigation on dibaryon structure made in 
the framework of chiral SU(3) model [7-9], in which the coupling between chiral fields 
and quarks is considered to describe the non-perturbative QCD effect, and therefore the 
huge data of nucleon-nucleon scattering phase shifts and hyperon-nucleon cross sections 
can be reproduced correctly. This investigation provided the di-omega (f2f2) + as the 
most interesting structure [8]. It is deeply bound with the binding energy about lOOMeV 
and its mean life time is in order of 10 -10 second because it can only undergo through 
weak decay. This charged dibaryon could be easily identified in experiments carrying on 
ultrarelativistic nucleus- nucleus collsions[8]. 

In order to find the deeply bound state (f2f2) + in RHIC experiments, we calculate 
its yield per event in central Au+Au collisions at y / J/v r /v — 130GeV. As we will see in 
the following, the hydrodynamic equation with transverse flow is used and measured data 
constitute inputs in all calculations to make results as reliable as possible. Cross sections 
for (f2f2) + yields in baryon-baryon scatterings are taken into account. 

2. Hadron distributions 

The hadron spectra measured in RHIC experiments allow us to analyse hadron 
distributions evolving with time via the Cooper- Frye formula [10], hydrodynamic equation 
and baryon number conservation. Assume hadron distributions are in the Boltzmann form 

f(X,p,T) = g\(r,r ] ,r)e^ T (1) 

where g, A, p and T are the degeneracy factor, fugacity, four-momentum and tempera- 
ture, respectively. The u is the four-velocity of fluid flow given by the cylindrical polar 
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coordinates r, and z [11,12]. The fugacity A depends on the proper time r, space- 
time rapidity r\ and r. The momentum distribution of particle given by means of the 
Cooper-Frye formula at freeze-out time Tf Q is 



mj_dm±dy (2ir) 2 Jo J»? min Jo 



m± COsh(y - r? ) e -[™xcosh(j / -„)-p ± ^cos0]/T v ^2 ^ 

where N, y, p±, m±_ and v r are respectively the number, rapidity, transverse momen- 
tum, transverse mass of hadron and transverse velocity of fluid flow. The -R(t/ d ) is the 
transverse radius of hadronic matter at freeze-out. The space-time rapidity spans from 
a minimum ?7 m i n to a maximum i] max . Any measured spectra of hadrons allow one to 
determine fugacities and transverse velocity at freeze-out. To gain the velocity profile at 
any proper time r, the hydrodynamic equation is used, 

dpT 1 * = (3) 

where the energy-momentum tensor is 

= (e + P)u^u v + Pg^ (4) 

The e and P are energy density and pressure, respectively. Because the dominant ingre- 
dient of hadronic matter is assumed to be pions, the equation of state is e = 3P. With 
the hydrodynamic equation for a boost-invariant longitudinal expansion along the z axis 
and a cylindrically symmetric transverse expansion [11], the transverse velocity v r (r,r) 
at any time prior to freeze-out is obtained. 

The dependence of hadron fugacity on time is obtained by the number conservation 

d^nu") = (5) 

which is followed by [13] 

dtX + Vrdr \ + 4^(7^ 3 ) + ^ r ( 7 T 3 ) + Xd r v r + A(^ + h = (6) 
7_z ^ 7 J J r t 



where 7 = 1/ yl~ 

The equation for fugacity has the divergent term Xv r /r as r — > 0, which causes 
difficulty in solving numerically the equation and leads to nonphysical solutions. There 
are four cases on curing this divergence: 

(1) v r ~ r n as r — > with n > 1; 

(2) A ~ r m as r — > with m > 1; 

(3) f r ~ r n and A ~ r m with n < 1 and m < 1 but n + m < 1 while r — > 0; 

(4) The term ^-(9 r (7T 3 ) cancels this term ^ as r — > 0, 

which has the analytic solution 
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T=(£)* (8) 



3 




The constant C\ is positive and in general depends on time. This solution becomes 
unphysical if it is singular as r — > 0. To avoid this, mandatorily force 7r = C<i with C2 
being nonzero and positive. The consequence is 

(9) 

Then the condition C 2 > r is required to give real transverse velocity. 

Before doing any calculations, we assume the hadronic matter freezing out at tj d = 
lOfm/c has a transverse radius R{rf ) = lOfm. The fugacity takes a form independent of 
77 and has the same dependence on r for all sorts of hadrons except for different normal- 
izations to fit rapidity densities. Given some trial functions for fugacities and velocities, 
the transverse mass spectra of 7r~, K~ and p at mid-rapidity in central Au+Au collisions 
at yj snn = 130GeV in Ref. [14] are reproduced. We find the following parametrizations 
showing r-dependences of fugacities and transverse velocity profiles at freeze-out 

\ 4 -o.7r ,9.65 + r 3 

A /o ocre u ', v r = { — - y (10) 

in the case (2) and 



e 0.2r 



in the case (4). Inserting the fugacities and transverse velocities into Eq. (1), we obtain 
hadron distributions as r — > 0: / — > 0. Evolve the hydrodynamic equation and the 
fugacity equation (6) back to a previous time from the freeze-out hypersurface, we still 
have / — > as r — > while the transverse velocity is kept as r — > as nonzero for the case 
(2) and the analytic function in Eq. (9) remains with a r-dependence of C2 for the case (4). 
Such / conflicts with the finite amplitude of a real hadron distribution as r — > 0. The cases 

(2) and (4) are therefore given up. Undoubtedly, nonzero hadron distributions beginning 
with the formation of hadronic matter may be found by seeking for a specific numerical 
method. But we do not consider here this task since we have no intention of finding all 
solutions of the hydrodynamic equation and the fugacity equation in coordinate space 
and all hadron distributions via Eq. (2) on the feeze-out hypersurface. Parametrizations 
for the fugacity and transverse velocity with the small r behavior described in the case 

(3) are not found. 

In the case (1) the fugacity and transverse velocity are not unique. The first set is 
, J 1, r < 7.5 r 

A ' o0C \ e-<-™> 2 , r>7.5 ^=104 (12) 

The v r depends linearly on r like Ref. [15]. The second and third sets of parametrizations 
are 

r 

v r = tanh (13) 
A /o oc I gV-T), r>l Vr = tanh 2l ( 14 ^ 
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In fact, the third set is obtained from the second set by making the replacement r — > 2r. 
For r > R{Tf ) the fugacity in the exponential form gives very small value and has 
negligible contribution to the integration over r. Then J R(Tfo) dr = ^ R(Tfo) dr is a very 
good approximation and the two sets give the same dependence of d 2 N/m±dm±dy on the 
transverse mass upon normalization. Actually, the substitution r — > ar with a > 1 leads 
to a class of fugacities and transverse velocities 

e V-T), r >£ v r = tanh — (15) 

which fits well the measured data of d 2 N/mj_dm±dy versus m± for ir~, K~ and p. Since 
the three curves in Fig. 1 of Ref. [14] are plotted in arbitrary units, the rapidity density 
<g = J™m ± dm ± m± d ^ ±dy is fitted to that generated in HIJING/BB [16] for p, p, n 
and n, respectively. In this way the magnitude of fugacity is completely fixed and the 
minimum and maximum of r] are individually equal to -4 and 4 for nucleon, -2.8 and 2.8 
for antinucleon. Meanwhile, a common freeze-out temperature T = 0.1 GeV is obtained 
for 7r~, K~ and p as extracted in Ref. [14]. 
Define an average velocity 

_ C T >°> drvr 

tv= ^r (16) 

The three sets of transverse velocity profiles in Eqs. (12)-(14) give v r w 0.48c, 0.67c, 0.83c, 
respectively. 

As already mentioned the measured data at RHIC energy might not well be re- 
produced in some theoretical models. The produced Q spectra is unknown. Even the 
knowledge of fi~ spectra measured in CERN-SPS experiments is very limited. To have a 
good sense on the fl~ spectra, we utilize the experimental data on the ratio of antiproton 
yield at RHIC to that at CERN-SPS. Since both the antiproton and Q~ at mid-rapidity 
can be created from junction and pomeron [16,17], we assume the ratio of Q~ yield at 
RHIC to that at CERN-SPS is the same as the ratio for antiproton. For rapidity densities 
of Q~ and antiproton, 

dN n - RHIC dN$i-,sPS _ dNp t RHic j dNp^ps 
dy dy dy dy 

It follows that 

dN n - jRHIC ri dNp^ RHIC , s 

-^— = c ^r- (18) 

where 

c = dN n ~ tSPS i dNp )S ps / ig x 
dy dy 

is taken as a constant and determined at mid-rapidity by 

c = r iy dN v ps i r *^ (20) 

i-0.5 dy 7-o.5 dy 

If the numerator is estimated from the values measured by WA97 Collaboration [18] and 
the denominator by NA49 Collaboration [19], C is about equal to 0.01. 
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3. Number of (Ntt) 2+ 



The chiral £77(3) quark model with the scalar nonet fields and pseudoscalar nonet 
fields has been presented to study (f2f2) + productions in baryon-baryon scatterings [20,21]. 
To create such a deeply bound state, the electromagnetic process 

fi + fi^(fifi) ++7 (21) 

and the mesonic process 

n + n -> (tttt) 0+ + 7] (22) 

are considered first. Since the state (f2f2) + is directly produced from two Q baryons, the 
two reactions are called direct production processes. In order to understand phyiscs well, 
the cross sections in Ref. [20] are replotted as functions of the center-of-mass energy of 
the two incoming omegas, y/s, in Fig. 1. 

The chiral SU(3) quark model shows a weakly bound state (iVT2) 2 + composed of 
a nucleon and an omega. The (NQ) 2 + dibaryon with spin 2 and positive parity can be 
formed through the electromagnetic process 

fi + iV ^ (M]) 2+ +7 (23) 

and the mesonic process 

Q + N -> (Ntt) 2+ + 7T (24) 

This (7Vfi) 2 + state enhances the yield of (QQ) + by three-quark exchange processes during 
its collisions with fl baryons 

Q + (Ntt) 2 + -> (tttt) + + N (25) 

Such (JlQ)o+ dibaryon is yielded from the collision of Q with nucleon and the accom- 
panying product (NQ) 2 +. The reactions are called indirect production processes. For 
present calculations it is better to redraw the cross sections in Ref. [21] as functions of 
the center-of-mass energy of omega and nucleon in Fig. 2 and of omega and (NQ) 2 + in 
Fig. 3. 

The number of (NQ) 2 + contained in hadronic matter at a proper time r is 
n (nq) 2+ = J d 4 x J ^-^f n (\ n ,p n ,T) J ^^f N (^N,PN,T)v rel a Nn (^) (26) 

which has the volume / d 4 x = J 27r d<p J^ f drr J^ T ^ drr /^™^ x dr] and where v reX is the relative 
velocity of two scattering particles and <j N si is the sum of cross sections for the reactions 
(23) and (24). The number of dibaryon (NQ) 2 + is calculated for a central Au+Au collision 
at y / J/viv — 130GeV. Results are plotted in Fig. 4. We take the hypothesis that (NQ) 2 + 
thermalizes immediately in hadronic matter once it is created. Then the momentum 
distribution of (iVf2) 2 + satisfies the Boltzmann form and the corresponding fugacity is 
determined by 

gj- f-R(~r) fflmax r^TT rOO r + OO 

N( N n) 2+ (r) = -j—-^ dr dr) d<p dm ± dyr\ N n(r,r],r) 

^Z7TJ JO Jrimin JO JmNSl J-oo 

ml cosh(y - r])e~ [m± cos Hy-v)-P±vr cos <t>]/T^fl r ^ 
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where m Nn is the dibaryon mass. The dibaryon fugacity Xnq takes the same r-dependence 
as a baryon except for normalization. 

4. Di-omega productions per event 

While the momentum distributions of Q and (NQ) 2 + are known, the di-omega yield 
per event is calculated via the formula 

N = J d4x J ^/i(Ai,Pi,T) J ^y 3 f2(X2,P2,T)v rel ann(V~s) (28) 

where / d 4 x = J 27r d(f) f^J° drr drr f%™^dr) with tj = 1.5fm/c corresponding to the 
temperature T m 0.19GeV. At the moment of forming hadronic matter, the interacting 
glod-gold system has already expanded. Presumably, the R(t) grows linearly with the 
proper time from R(tj) = 7fm to -R(r/ ) = lOfm. The ann include cross sections for the 
reactions (21), (22) and (25). 

Given a fugacity and transverse flow velocity, number of dibaryon (QQ) + yielded 
in a central Au+Au collision at y^s^N — 130GeV is listed in Table 1. N 1 , N 2 and A 5 are 
individually the number of (f2fi) + created through the reactions Q + Q — > (fifi) + + 7, 
f2 + f2 — > (f2fi) + + i] and f2 + (Af2) 2 + — > (f2fi) + + N. The total number obtained in the 
work is 

N tot = N 1 + N 2 + N 5 (29) 
where Ai, N 2 and A 5 are calculated with the use of Eq. (28). 

5. Discussions and conclusions 

We have shown that different hadron distributions at freeze-out lead to different 
(f2f2)o+ numbers produced in a Au+Au collision and have obtained that 7r~, K~ and p 
freeze out at the same temperature T = O.lGeV while the radius of transverse extension of 
cylindrically symmetric hadronic matter is lOfm. If the freeze-out temperature or trans- 
verse radius changes, these distributions vary. The fugacities and transverse velocities in 
Eqs. (14) and (15) result from the replacement r — > 2r and r — > ar in the parametrization 
(13). However, the hydrodynamic equation and the fugacity equation (6) are not invariant 
under the replacement. This is one reason why the number of (£7f2)o+ corresponding to 
the parametrization (13) is different from that to the other parametrization (14). Another 
reason is that a factor a 2 accounts for the variation of the integration / drr under the 
replacement r — > ar. 

In the indirect production processes the number of (NQ) 2 + depends on the freeze- 
out temperature, the final transverse radius and the final proper time when the hadronic 
matter becomes free particles. The yield of (f2f2) + in the indirect production processes 
is actually proportional to the square of the distributions of proton and neutron. Fig. 4 
exhibits very small number of (NQ) 2 +. However, the cross sections for each of indirect 
production processes, especially the exothermal reaction Q + (NQ) 2 + — > (f2fi) + + N in 
Fig. 3, are much larger than the cross sections for the direct production processes. The 
amount of di-omega gained in the indirect production processes eventually exceeds that 
in the direct production processes. Nevertheless, this is true only near freeze-out when 
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the (NQ) 2 + number is relatively larger, and false near hadronization when the (iVT2) 2 + 
number is lack. 

It is shown by the three curves in Fig. 4 that the (NQ) 2 + dibaryon number increases 
rapidly during the time interval from r = 1.5fm/c to 3fm/c. Almost all of (NQ) 2 + are 
produced in the early time of hadronic matter and sensitive to early-time thermalization. 
The variation of dibaryon number with time is controlled by high baryon densities in the 
early time of hadronic matter and later on low baryon densities. We are also interested 
in calculating the dependence of (fifi) + number on the proper time with, for instance, 
the fugacity and transverse velocity in Eq. (12). Unlike the (NQ) 2 +, we find that the 
di-omega number changes slowly with time and almost all of (f2f2) + are produced in the 
late time of hadronic matter. Even if the (NQ) 2 + number increases rapidly, its density 
rises up not fast because the matter volume increases. Together with the reducing density 
of £l~ with time, the (fif2) + is produced slowly in hadronic matter. 

The reaction Q + Q — > (fifi) + + r\ has the threshold energy 3.776GeV higher than 
the value 3.344GeV for il + VI — > (f2f2) + + 7. Peaks of cross sections are shown in Fig. 
1 for both reactions. Thereby the f2 has larger average energy in the rj yield process 
than in the 7 yield process. This considerably reduces the (f2fi) + yield from the former 
reaction compared to the latter by the momentum distribution of Q. However, the ratio 
N 2 /N\ m 2.5 in Table 1 reverse the case. This comes out of the following factors: the 
reaction for rj production has larger cross section and a wider peak than for 7; higher 
relative velocity of two f2s is demanded for 77 production; d 3 pid 3 p 2 in Eq. (28) is greatly 
enhanced by the larger average momentum of Q needed for i] yield. Therefore, more 
(f2f2) + dibaryons are produced in the rj yield process than in the 7 yield process. 

We have not estimated di-omega yields in the expansion process from the moment of 
initial nucleus-nucleus collision to the time 7/ = 1.5fm/c when hadronic matter appears. 
During this period the central Au+Au collision at y/s NN = 130GeV may lead to decon- 
fined gluons and quarks in nonequilibrium, if we believe at CERN-SPS energy the Pb+Pb 
collisions are on the boundary of phase transition involving hadronic and partonic degrees 
of freedom. While the quarks and gluons hadronize at 77, it is possible for six strange 
quarks to combine into a dibaryon (Qfi) + as well as six light quarks to form (NQ) 2 +. 
Then the number of (NQ) 2 + at 7/ = 1.5fm/c is not zero and moves upwards these curves 
in Fig. 4. 

It is unavoidable to address in-medium effect on the production of (f2f2) +. The 
mass of strange quark is reduced and the quark-quark potential is screened in hadronic 
matter at finite temperature. Obviously the screening decreases the binding energy of 
(f2f2) +. However, the kinetic-energy contribution increases the binding energy by the 
dropping mass of strange quark. We can still expect an appreciable binding energy for 
(f2f2) + in hadronic matter with reference to the large binding energy of free (f2f2) +. In 
addition, hadronic matter modifies the relative-motion wave function of two omegas inside 
the dibaryon and the coupling of rj to Q in Q + Q — > (f2f2) + + 77 and the coupling of rr 
to N in f2 + N — > (iVT2) 2 + + n. These result in different cross sections for the direct and 
indirect production processes in comparison to Figs. 1-3 and different (fifi) + number 
yielded relative to Table 1. Exploring in-medium effect remains in a future work. 

In conclusion, we have obtained the number of dibaryon (f2f2) + produced in a central 
Au+Au collision at ^s NN = 130GeV is of the order of 1CT 9 ~ 1CT 7 . Even if the number 
is very small, enough events will be accumulated to observe the yield in the coming years 
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of running of RHIC machine. Besides, we have shown: the mesonic processes have larger 
contributions than the electromagnetic processes; at feeze-out the indirect production 
processes have larger contributions than the direct production processes; di-omegas are 
mostly produced in late-time hadronic matter. 
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Figure captions 



Fig. 1. The left and right panels show cross sections for Q + Q — > (f2f2) + + 7 and 
Q + Q — > (f2f2) + + 77, respectively, as functions of the center-of-mass energy of two 
incident omegas. 

Fig. 2. Dotted, dot-dashed and solid lines in the left panel show cross sections for 
Q + p — > (jof2) 2 + +7, CI + n — > (nfi) 2 + + 7 and their sum, respectively. The right panel 
shows the cross section for f2 + N — > (A^fi) 2 + + 7r where iV includes proton and neutron. 

Fig. 3. Cross section for fl + (NQ) 2 + — > (fir2) + +iVasa function of the center-of-mass 
energy y/s. The includes proton and neutron. 

Fig. 4. Solid, dashed and dot-dashed lines are the numbers of (NQ) 2 + produced with the 
parametrizations for fugacities and transverse velocities in Eqs. (12)-(14), respectively. 
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Table 1: Number of (f2f2) + produced per event with the three sets of fugacities and 
transverse velocities in Eqs. (12)-(14). 
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